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QUANTUM  EXTENSION  OF  CHILD-LANGMUIR  LAW 


The  Child-Langmuir  lav^  gives  the  maximum  electron  current  that  can  be 

transmitted  across  a  parallel,  planar  gap  in  terms  of  the  incident  energy 

of  the  electrons  and  the  gap  bias  voltage.  This  maximum  value,  known  as 

the  limiting  current,  arises  because  the  space  charge  in  the  diode  presents 

a  potential  barrier  to  the  incident  electrons.  While  there  are 

modifications  due  to  geometrical  and  relativistic  effects,  the  limiting 

current  remains  a  fundamental  quantity  characterizing  the  beam-gap 
.  .  2 

interaction. 


In  the  emerging  fields  of  nan-.-electronics,  tunneling  microscopy  and 
vacuum  micro-alectronics,  diode  gaps  and  junctions  with  scales  down  to  tens 
of  angstroms  are  being  considered.  On  such  a  microscopic  scale,  quantum 


effects  may  no  longer  be  neglected.  Several  questions  then  arise:  Is 
there  a  limit  on  the  current  that  can  be  transmitted  across  the  gap  when 


quantum  effects  are  taken  into  account?  Is  the  transmitted  current 

quantized?  How  is  the  classical  value  recovered,  in  conformity  with  the 

correspondence  principle?  This  paper  addresses  these  questions. 

Here,  we  extend  the  classical  work  of  Child  and  Langmuir  to  the 

quantum  regime  by  considering  a  parallel,  planar  gap.  We  will  use  the 

familiar  mean-field  theory  expressed  by  the  self-consistent,  coupled 

Schrodinger  and  Poisson  Equations  in  the  Hartree  approximation.  Similar 

approaches  have  been  taken,  for  example,  to  study  the  effects  of  space 

3 

charge  on  the  device  characteristics  of  superlattice  structures.  Thus,  we 


□  □ 


and  the  Poisson  Equation 


«4  -  .♦♦*/«.  (2) 

dxz  0 

in  the  gap  region,  defined  to  be  0  <  x  <  D. 

Implicit  in  Eqs.  (1)  and  (2)  are  the  following  assumptions,  (a)  The 

usual  Schrodinger  wave  function  <p(x)  is  interpreted  as  the  density  of  a 

continuum  electron  fluid,  Llit  .:uieroscopic  numbet  density  being  u  =  • 

(b)  The  electrostatic  potential  V(x)  is  time-independent,  thus  allowing  the 

separation  of  the  temporal  dependence  in  the  form  ^(xje-*^*^1.  (c)  All 

electrons  are  non-relativistic,  and  they  enter  the  gap  at  x  =  0  with  the 

same  kinetic  energy,  (d)  The  electron  density  in  the  gap  is  sufficiently 

small  that  single  particle  wave  functions  do  not  overlap  significantly. 

(e)  Any  one  electron  sees  an  average  electric  field  due  to  all  others 

present  in  the  gap.  Thus,  vhat  we  are  considering  is  a  situation  where  the 

electron  density  in  the  gap  is  sufficiently  high  that  the  self  fields  of 

the  electrons  are  no  longer  negligible,  but  is  low  enough  that  we  may  omit 

consideration  of  the  exclusion  principle.  If  the  average  number  of 

electrons  becomes  too  low,  large  fluctuations  in  both  the  current  and  the 

self-potential  will  occur,  and  the  quantum  correction  calculated  in  this 

paper  should  then  be  viewed  only  as  a  first  approximation. 

The  boundary  conditions  to  Eqs.  (1)  and  (2)  are  formulated  as  follows. 

Let  V  be  the  gap  bias  voltage,  which  may  be  zero,  positive  or  negative. 

S 

Then  V(0)  «=  0  and  V(D)  =  V  .  The  boundary  conditions  on  the  wave  function 

6 

are  derived  from  the  conditions  that  «|>  be  matched,  at  x  =  0,  to  the  sum 

ik.x 

of  an  incident  plane  wave  <J^(x)  =  Ae  and  a  reflected  plane  wave  «|»r(x)  = 

-ikjX  ^2X 

Be  and,  at  x  =  D,  to  a  transmitted  plane  wave  ^(x)  =  Ce  .  Here, 


2 


A,  B,  C  are  constants  and  k^  =  (2mE/fi2)^2,  k2  =  [2m(E+eVg)/fi2] *^2.  Charge 

conservation  requires  that  the  current  density  J  =  e(if»/2m)(<H'*’  -  )  be 

constant  for  all  x.  Here,  a  prime  denotes  a  derivative  with  respect  to  x. 

In  terms  of  the  wave  amplitudes,  J  =  (efi/m^  |C |2  =  (efi/nOk^  |A  |2  -  |b|2). 

In  keeping  with  the  classical  theory,  we  shall  assume  that  E  and  V  are 

& 

given,  with  E  +  eV  >  0.  We  shall  determine  the  conditions  on  J  for  the 

o 

existence  of  solutions  to  Eqs.  (1)  and  (2). 

To  see  clearly  the  transition  from  the  classical  to  the  quantum 

regime,  we  find  it  convenient  to  use  non-dimensional  quantities: 

x  =  x/D,  V  =  eV/E,  E  =  E/eV  .  J  =  J/J  ,  n  =  n/n  =  |<J»|2/n  . 

s  s  s  s 

2  2 

♦  =  eV  /E,  where  the  voltage  scale  V  =  “f:  /(2emlJ  ) ,  the  current  density 

o  o  S 

scale  Jg  =  eo'ti"^/t4m2eD^),  and  the  number  density  scale  ng  «=  eQ^2/(2e2mD^) . 

__  I/O  __  iG(x^ 

We  next  represent  the  wave  function  y(x)  =  (n  E)  p(x)e  v  '  in  terms  of 

s 

the  non-dimensional  amplitude  p(x)  and  phase  0(x),  both  assumed  real. 
Equations  (1)  and  (2)  yield  the  following  coupled  equations  for  p(x)  and 
V(x): 


1  d2p 

=  “Zf  + 

E  dxz 


(1+  V)  - 

P 


P  =  0 


2 

P 


(3) 

(4) 


where  we  have  introduced  the  dimensionless  "perveance" 

A  *  2J/E3/2  (5) 

which  is  proportional  to  the  current.  In  terms  of  A,  the  phase  0(x)  is 
given  by 


3 


e(x) 


X 

[  dx/p2(x) 


(6) 


where,  without  loss  of  generality,  we  have  assigned  0(1)  =  0  [Ref.  4].  The 

4 

boundary  conditions  to  Eqs.  (3)  and  (4)  are 


V<0)  =  0  (7) 

V(l)  =  4>g  (6) 

p(l)  =  (X/4)1/2/(l  +  *g)1M  (9) 


P’(l)  =  0.  (10) 

It  can  easily  be  shown  that  the  classical  limit  is  obtained  by  simply 

n  __ 

ignoring  the  first  term  P  /E  in  Eq.  (3)  and  that  the  quantum  behavior 
enters  only  through  that  term.  The  classical  behavior  dominates  when 
E  »  1,  and  the  transition  to  the  quantum  regime  is  expected  to  occur  when 
E  =  0(1).  Thus,  in  the  classical  limit  (E  -*  “},  Eq.  (3)  gives 


(1  +  'r* 


which,  when  cast  back  in  dimensional  form,  is  simply  the  statement  of 
energy  conservation  in  the  classical  description  of  electron  motion. 
Indeed,  the  Child-Langmuir  law  may  be  recovered  by  substituting  Eq.  (11) 
into  Eq.  (4).  The  resultant  second  order  differential  equation  in  V, 
subject  to  boundary  conditions  (7)  and  (8),  may  be  shown  to  admit  no 
solution  whenever  X  >  Xc>  where 

xc  =  r  t1  +  (1  +  V1/2]3-  (12) 


4 


1  2 

This  is  essentially  the  Child-Langmuir  law  ’  in  normalized  form.  It  gives 
the  maximum  current  that  can  be  transmitted,  in  steady  state,  in  terms  of 
the  injection  energy  E  and  the  gap  bias  voltage  V  .  Note  that  X£  is 
independent  of  E  [Fig.  1]. 

Retaining  the  p"  term  in  Eq.  (3),  our  problem  becomes:  for  specified  E 
and  f  ,  for  what  values  of  X  do  equations  (3)  and  (4)  have  solutions, 

o 

4 

subject  to  Eqs.  (7 )— ( 10) .  We  have  obtained  numerical  solutions  for  a  wide 

range  of  values  of  E  and  ♦  .  In  particular,  we  have  found  that  there  is  a 

S 

critical  value  of  X,  called  X^,  above  which  no  solution  exists.  According 
to  this  formulation,  Xq  represents  the  maximum  current  that  can  reach  x  = 

1,  independent  of  the  nature  of  the  emitter  at  x  ■  0,  since  the  boundary 
conditions  (9)  and  (10)  specify  only  the  transmitted  flux  at  x  =  1. 

Our  results  are  summarized  in  Fig.  1,  which  shows  X^  as  a  function  of 
E  for  different  values  of  4  .  The  classical  value,  from  Eq.  (12),  is  also 

o 

shown.  Note  that  the  classical  value  is  indeed  approached  for  large  values 

of  E  vhile,  for  small  values  of  E,  X^  greatly  exceeds  the  classical  value. 

We  attribute  this  finding  to  the  tunneling  of  electrons  through  the 

potential  barrier  presented  by  the  average  space  charge  field  of  other 

particles  in  the  gap.  In  fact,  using  the  numerically  computed  potential 

barrier  (-V),  [Fig.  2],  we  have  found  that  the  WKBJ  estimates  on  the 

tunneling  across  such  a  barrier  are  consistent  with  the  transmission 

coefficient  that  was  computed  numerically  [Fig.  3].  For  small  E,  Fig.  1 

suggests  X  «  1/E. 

9 

The  transmittable  current  is  not  quantized  in  the  above  formulation. 
For  all  X  <  X^,  solutions  to  Eqs.  (3),  (4)  subject  to  (7)  -  (10)  can  always 
be  found. ^  Figure  2  shows  a  sample  solution.  From  the  numerical 


5 


2 

solutions,  we  may  construct  a  reflection  coefficient  CR  *  l^r^i |  and  a 
transmission  coefficient  C^,  m  (k2/k^)  1^/^  |  relative  to  the  incident  flux 


(current) . 
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Clearly,  CR  +  CT 


1.  The  transmission  coefficient  is  shown 


in  Fig.  3. 

0  1 6  _3 

As  an  example,  take  D  =  30A,  then  V  =  4.24mV,  n  =  2.6  x  10  cm  ,  J 

s  s  s 

2  — 

=  8.06  kA/cm  .  If  we  further  take  E  =  1  and  ♦  =  0.5,  then  Fig.  1  gives  X 

&  Q 

=  102,  whereas  the  classical  theory  Eq.  (12)  gives  X£  =  19.6.  This  value 


of  Xq  means  that  the  maximum  current  density  that  can  be  transmitted  across 

5  2 

such  a  gap  is  4.11  x  10  A/ cm  from  the  quantum  mechanical  theory.  On  the 


other  hand,  according  to  the  classical  theory,  the  maximum  current  density 

4  2 

would  only  be  7.89  x  10  A/cm  ,  a  factor  of  5  lower.  The  solutions  for  this 

example  with  X  =  95  are  shown  in  Fig.  2,  from  which  we  can  obtain  the 

average  electron  density  <n>  and  the  spatial  scale  L  over  which  the  wave 

1/3 

function  varies.  One  can  readily  deduce  that  <n>  L  <  1.  Thus,  the 
electron  density  is  sufficiently  low  to  ignore  the  exclusion  principle, 
justifying  the  present  mean  field  theory  approach. 

In  summary,  self  consistent  solutions  to  Eqs.  (1),  (2)  have  been 
constructed.  From  these  solutions,  we  have  found  that  the  macroscopic 
current  that  can  be  transmitted  across  a  gap  can  exceed  the  classical  value 
of  Child  and  Langmuir,  sometimes  by  a  large  amount,  because  of  tunneling 
effects.  Dimensionless  parameters  have  been  identified  whereby  the 
transition  from  the  quantum  regime  to  the  classical  regime  may  be  expected. 
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on  the  gap,  whereas  (9)  and  (10)  follow  from  the  requirement  that  the 
solution  at  x  ■  1  matches  a  (preassigned)  transmitted  wave  <|<t(x).  To 
integrate  (3)  and  (4),  we  use  (8)  -  (10)  and,  in  addition,  assume  a 
value  for  V'(l),  as  initial  conditions  at  x  =  1.  The  value  V'(l)  is 
adjusted  sc  that  condition  (7)  is  satisfied,  after  integrating  (3), 

(4)  back  to  x  =  0.  The  incident  wave  and  the  reflected  wave 
amplitudes  can  then  be  inferred  from  the  numerical  solutions. 

5.  For  A  <  \q,  we  find  two  solutions  of  Eqs.  (3),  (4)  satisfying  (7)  - 

(10),  for  specified  values  of  E,  ♦  ,  and  X.  For  various  reasons,  we 

argue  that  the  one  with  higher  potential  energy  is  inaccessible.  In 
this  paper,  we  focus  only  on  the  solutions  with  lower  potential 
energy.  When  X  =  X^,  the  two  solutions  merge. 
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6. 


The  present  paper  essentially  treats  a  nonlinear  scattering  problem  of 
beam  injection  into  a  gap  by  an  external  source.  It  is  conceivable 
that  a  different  physical  situation  would  require  different  boundary 
conditions  on  ♦  at  x  =  0  that  would  lead  to  the  interesting 
possibility  of  quantization  of  both  E  and  X.  The  existence  of  such 
states,  and  their  stability,  will  be  the  subjects  of  a  future 
publication.  We  should  add,  however,  that  regardless  of  the  boundary 
conditions  on  V  that  would  be  imposed  at  x  =  0,  the  normalized 
critical  current  X^  calculated  in  this  paper  is  still  the  upper  limit. 
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